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Abstract 

We report the Fermi Large Area Telescope (LAT) detection of a 7 -ray source at the position of SAX 
J1808.4—3658. This transient low-mass X-ray binary contains an accreting millisecond pular, which is 
only seen during its month-long outbursts and likely switches to be rotation powered during its quiescent 
state. Emission from the 7 -ray source can be described by a power law with an exponential cutoff, the 
characteristic form for pulsar emission. Folding the source’s 2.0-300 GeV photons at the binary orbital 
period, a weak modulation is seen (with an H-test value of ^17). In addition, three sets of archival XMM- 
Newton data for the source field are analyzed, and we find only one X-ray source with 3-4tT flux variations 
in the 2tT error circle of the 7 -ray source. However based on the X-ray properties, this X-ray source is 
not likely a background AGN, the major class of Fermi sources detected by LAT. These results support 
the possible association between the 7 -ray source and SAX J1808.4—3658 and thus the scenario that the 
millisecond pulsar is rotation powered in the quiescent state. Gonsidering a source distance of 3.5 kpc for 
SAX J1808.4—3658, the 0.1-300GeV luminosity is 5.7x10^^ergs“^, implying a 7 -ray conversion efficiency 
of 63% for the pulsar in this binary. 
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1. INTRODUCTION 

Millisecond radio pulsars (MSPs) are formed from neu¬ 
tron star low-mass X-ray binaries (LMXBs; Bhattacharya 
& van den Heuvel 1991). The primary neutron star in an 
LMXB can gain sufficient angular momentum by accreting 
material from the companion through an accretion disk, 
and thus be ‘recycled’ to a spin period of milliseconds. 
The discovery of millisecond X-ray pulsations in the tran¬ 
sient LMXB SAX J1808.4—3658 has confirmed the for¬ 
mation scenario from the observational side (Wijnands & 
van der Klis 1998; Chakrabarty & Morgan 1998). Thus far 
over a dozen of so-called accreting millisecond X-ray pul¬ 
sars (AMXPs) have been found (Patruno & Watts 2012). 
Nearly all of them are in transient systems, and for these 
transients. X-ray pulsations are seen only during their X- 
ray outbursts. 

Interestingly, it was suggested that several AMXPs ac¬ 
tually switch to be rotation powered pulsars during their 
quiescent states (Burderi et al. 2003; Wang et al. 2013 
and references therein), although no direct evidence, such 
as pulsed radio emission (Burgay et al. 2003), was found 
over the time. The recent observational identification of 
the MSP binary J1824—24521 in the globular cluster M28 
has firmly confirmed the suggestion (Papitto et al. 2013). 
The binary was observed to have an X-ray outburst, and 
during the outburst, the previously known radio MSP in 
the binary switched to appear like a typical AMXP. This 
confirmation has thus identified an interesting feature for 
the evolution from LMXBs to MSP binaries, and we may 


suspect that either these systems would probably be at 
the end of their LMXB phase or it could be a common 
feature for the quiescent state of transient neutron star 
LMXBs (e.g., Heinke et al. 2014). 

The similar type of feature has also been seen 
in the recently identihed two transitional MSP bi¬ 
naries: J1023-I-0038 (Archibald et al. 2009) and 

XSS J12270—4859 (Bassa et al. 2014 and references 
therein). Extensive observational studies have shown that 
they can switch between the states of having an accretion 
disk and being disk free. One particular property in them 
is that they have sufficiently bright 7 -ray emission and 
are detectable by the Fermi Large Area Telescope (LAT; 
see Tam et al. 2010; Stappers et al. 2014; Takata et al. 
2014 for J1023-I-0038; and see Hill et al. 2011; de Martino 
et al. 2013; Xing & Wang 2014 for XSS J12270-4859). 
The emission is variable, stronger in the active accre¬ 
tion state than that in the disk-free state (Stappers et al. 
2014; Takata et al. 2014; Xing & Wang 2014). Given the 
property similarities between the two MSP binaries and 
AXMP binaries, it is thus highly possible that AMXPs 
would also have significant 7 -ray emission, as they stay in 
their quiescent state and would be rotation-powered most 
time. This possibility has been explored by Xing & Wang 
(2013) by searching for 7 -ray emission from four AMXP 
systems, which include SAX J1808.4—3658. Nearly four 
year LAT data for the four AMXPs were analyzed, but no 
7 -ray emission was found. However, given the improved 
sensitivity of Fermi over the last two years (see § 3.1 for 
details), we re-analyzed the LAT data for them. We found 
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significant 7 -ray emission at the position consistent with 
that of SAX J1808.4—3658. Here in this paper, we report 
the results. 

2. Observation 

LAT is a 7 -ray imaging instrument onboard the Fermi 
Gamma-ray Space Telescope. It makes all-sky survey 
in an energy range from 20 MeV to 300 GeV (Atwood 
et al. 2009). In our analysis, we selected LAT events 
from the Fermi Pass 7 Reprocessed (P7REP) database 
inside a 20 ° x 20 ° region centered at the optical position 
of SAX J1808.4-3658, which is R.A.=18‘’08“27i62, 
Decl.=—36°58'4153 (equinox J2000.0) obtained in 
Hartman et al. (2008). We kept events during the time 
period from 2008-08-04 15:43:36 (UTG) to 2014-11-10 
21:04:57 (UTG) and in the energy range of 100 MeV to 
300 GeV. In addition, only events with zenith angle less 
than 100 deg and during good time intervals were selected. 
The former prevents the Earth’s limb contamination, and 
for the latter, the quality of the data was not affected by 
the spacecraft events. 

3. Data Analysis and Results 

3.1. Source Identification 

We included all sources within 20 deg in the Fermi sec¬ 
ond source catalog (Nolan et al. 2012) centered at the 
position of SAX J1808.4—3658 (Hartman et al. 2008) to 
make the source model. The spectral function forms of the 
sources are provided in the catalog. The spectral param¬ 
eters of the sources within 5 deg from SAX J1808.4—3658 
were set free, and all other parameters of the sources 
were fixed at their catalog values. A point source at 
the optical position of SAX J1808.4—3658 was also in¬ 
cluded in the source model, with its emission modeled 
by a simple power law. In addition, we used the spec¬ 
trum model glLiem_v05_revl.fits and the spectrum file 
iso_source_v05.txt for the Galactic and the extragalactic 
diffuse emission, respectively, in the source model. The 
normalizations of the diffuse components were set as free 
parameters. 

Using the LAT science tools software package v9r33p0, 
we performed standard binned likelihood analysis to the 
LAT data. Events below 200 MeV were rejected be¬ 
cause of the relative large uncertainties of the instrument 
response function of the LAT in the low energy range. 
Energy ranges of 0.2-300, 0.5-300, 1-300, and 2-300 GeV 
were tested in the analysis. A source at the optical po¬ 
sition was detected with Test Statistic (TS) values of 32, 
34, 26, and 31, respectively. The TS value at a specific po¬ 
sition, calculated from TS= —21og(Lo/Li) (where Lq and 
Li are the maximum likelihood values for a model without 
and with an additional source respectively), is a measure¬ 
ment of the fit improvement for including the source, and 
is approximately the square of the detection significance 
of the source (Abdo et al. 2010). Thus the source was best 
detected in 0.5-300 GeV with a significance of ~5cr. We 
extracted the TS maps of a 2° x 2° region centered at the 
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position of SAX J1808.4—3658 in the four energy ranges, 
with all sources in the source model considered except the 
source we found. No catalog sources are within the square 
region. In Figure 1, the 0.5-300 GeV TS map is shown. 

We ran gtfindsrc in the LAT software package to deter¬ 
mine the position of the source using photons in 0.5-300 
GeV, and obtained the best-fit position R.A. = 272?15, 

Deck = —37P06 (equinox J2000.0), with Icr nominal un¬ 
certainty of 0.°05. The 2cr error circle is marked in Figure 1 
as a dark dashed circle. The optical position of SAX 
J1808.4—3658 (mark by a dark cross in Figure 1) is 0.°08 
from the best-fit position and within the 2 cr error circle, 
suggesting possible association of the 7 -ray source with 
SAX J1808.4—3658. Below we considered the source as 
the candidate 7 -ray counterpart to SAX J1808.4—3658. 

In our TS maps, separate excess 7 -ray emission at the 
northwest corner appears (Figure 1). We investigated 
whether it possibly contaminated our detection of the can¬ 
didate counterpart. We found that it is consistent with be¬ 
ing a point source at a position of R.A. = 271.°4, Deck = 

—36?4 (equinox J2000.0; la nominal uncertainty is 0?1). 

Including this source in our source model, it can be totally 
removed from the TS maps, and the results of the position 
and spectrum (see Section 3.2) of the counterpart source 
did not have significantly differences (consistent within 
uncertainties). 

In addition, since the source position is toward the 
Galactic center direction (Gf, ~ — 8?1), we also checked 
if the uncertainty on the Galactic diffuse emission could 
produce false detection of the 7 -ray source. We manually 
increased the normalization of the Galactic diffuse com¬ 
ponent to a value 5cr above the best-fit value, the >0.5 
GeV 7 -ray emission at the position of SAX J1808.4—3658 
was still detected, with TS~32. We found only when we 
increased the Galactic diffuse emission by 10% (approxi¬ 
mately 77a above the best-fit value), the >0.5 GeV emis¬ 
sion was then detected with TS~9 (i.e., ^3a detection 
significance). 

The 7 -ray source was not detected in the previous 
search using nearly four-year LAT data in >0.2 GeV en¬ 
ergy range (Xing & Wang 2013). We re-analyzed the LAT 
data in the same time interval from 2008-08-04 15:43:36 
(UTG) to 2012 July 8 18:59:57 (UTG), and obtained a 
TS value of ~25 at the position of SAX J1808.4—3658, 
which is much higher than the value of ^^3 previously 
reported. Therefore the detection is due to the im¬ 
proved sensitivity of the Fermi telescope. The database 
used in Xing & Wang (2013) is Fermi Pass 7, compar¬ 
ing to Pass 7 Reprocessed in this work. The LAT sci¬ 
ence tools software package and the Instrument Response 
Functions (IRFs) have also been updated from v9r27pl 
to the current v9r33p0 and from P7SOURCE_V6 to 
P7REP_SOURCE_V15, respectively. In addition, the dif¬ 
fuse emission models have also been updated. All these 
changes have improved the point source detection sensi¬ 
tivity of the Fermi/LAT^. 

^ http: / / www.slac.stanford.edu/exp / glast / groups / canda/lat_Performance.htm 
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3.2. Spectral Analysis 

Including the 7 -ray source in the source model, we per¬ 
formed standard binned likelihood analysis to the LAT 
data, with emission of this source modeled with an ex¬ 
ponentially cutoff power law, which is characteristic of 
pulsars (Abdo et al. 2013). In addition a simple power 
law was also used. We used data in >0.2 GeV energy 
range to obtain a overall description of the 7 -ray spec¬ 
trum of the source. A photon index of P = 2.2±0.1 
with a TSp/ value of ~32 was obtained for the power- 
law model, and a photon index of P = 1.6±0.4 and a 
cutoff energy of = 5.5±3.7 GeV with a TSe^jp value 
of ~37 were obtained for the exponentially cutoff power- 
law model. The low energy cutoff was thus detected 
with > 2(7 significance (estimated from y^TScutoff, where 
TSouto// — TSexp — TSp/ ~ 5; Abdo et al. 2013). While 
the significance is low, this result also favors the possible 
association of the 7 -ray source with SAX J1808.4—3658. 
These spectral results are summarized in Table 1. 

We then extracted the 7 -ray spectrum for the source, by 
considering the emission as a point source with a power- 
law spectrum at the optical position of SAX J1808.4—3658 
and performing maximum likelihood analysis to the LAT 
data in 10 evenly divided energy bands in logarithm from 
0.1-300 GeV. The photon index was fixed at 2.2. Only 
spectral points with TS>4 were kept, and the 95% upper 
limits in other energy bins were derived. The spectrum 
extracted by this method is less model dependent and 
provides a more detailed description for the 7 -ray emission 
of the source. The obtained spectrum is shown in Figure 2, 
and the energy flux values are given in Table 2. It can be 
seen that the exponentially cutoff power law fits the data 
better, particularly in low energy ranges where no 7 -ray 
emission was significantly detected. 

3.3. Variability Analysis 

We performed timing analysis to the LAT data of the 
SAX J1808.4—3658 region to search for possible 7 -ray 
pulsations. We folded the LAT data according to the 
X-ray ephemeris given in Hartman et al. (2009). The 
optical position of SAX J1808.4—3658 was used for the 
barycentric corrections to photon arrival times, and pho¬ 
tons within Rmax {Rmax ranges from 0 .° 1 - 1.°0 with a step 
of 0.°1) from the position were collected. Different energy 
ranges (>0.2, >0.5, >1, and >2 GeV) were tested in fold¬ 
ing. No pulsation signals were detected, and the H-test 
values were smaller than 9 (corresponding to <3tT detec¬ 
tion significance; de Jager et al. 1989). 

We folded the LAT data using the orbital parameters 
given in Hartman et al. (2009). We found that the highest 
orbital signal was revealed in the >2 GeV energy range 
using photons within 0.°6 from the optical position of SAX 
J1808.4—3658. The folded light curve, which has an H- 
test value of ~17 (corresponding to > 3(7 detection signif¬ 
icance, de Jager et al. 1989), is shown in Figure 3. The 
phase zero is set at the ascending node of the pulsar in 
SAX J1808.4-3658. 

We also obtained the light curves for the 7 -ray source. 


with different time intervals (e.g., 30, 100, and 300 days) 
used. Due to the faintness of the source, no significant 
flux variations can be determined from the light curves. 

3 . 4 . XMM-Newton data Analysis 

We searched in the SIMBAD Astronomical Database. 
Most sources identified in the 7 -ray source region are not 
high-energy but star-type objects. There are two other X- 
ray sources, SAX J1808.5-3703 and SAX J1809.0-3659, 
previously reported in the region (Wijnands et al. 2002; 
Campana et al. 2002). While the two sources are 0.°02 and 
0?11 away from our best-fit position, respectively, there 
are also other X-ray sources detected in the region but 
not studied in detail (Gampana et al. 2002). 

We thus searched and found three archival XMM- 
Newton observations of SAX J1808.4—3658 available. 
They were carried out on 2001 Mar. 24 (ObsID ; 
0064940101; exposure 39.5 ks), 2006 Sept. 15 (ObsID 
: 0400230401; exposure 55.1 ks), and 2007 Mar. 10 
(ObsID ; 0400230501; exposure 57.8 ks). We analyzed the 
European Photon Imaging Gamera (EPIG) pn and MOS 
data using the standard tools of the XMM-Newton Science 
Analysis Software (SAS, version 14.0). For the first ob¬ 
servation, pn was used in the timing mode and the data 
were not included in our analysis. We excluded the high 
particle flaring background by creating the good time in¬ 
tervals (GTI) based on the count rate cut-off criteria. We 
extracted the full-held background light curve in the 10 - 
12 keV band and selected the GTI with count rate < 0.8 
and < 0.4ctss“^ for pn and MOS data respectively. The 
data were then hltered to the good X-ray events (FLAG 
== 0) with PATTERN < 4 for pn and PATTERN < 12 
for MOS in the 0.3-10 keV energy band. 

We then performed the source detection routine 
(EDETECT_CHAIN) on 2007 EPIG-pn data and identihed 
17 held X-ray sources in the 2cr Fermi error circle (ra¬ 
dius of 0?1) other than the AMXP, SAXJ 1808.4—3658. 
Their positions were astrometrically calibrated by corre¬ 
lating the sources detected in the whole pn held with 
the USNO Bl.O optical catalog (Monet et al. 2003; 
the SAS task EPOSCORR was used). The X-ray held 
is shown in Figure 4. We obtained the source counts 
of the 17 sources in both pn and MOS data of the 
three observations. They were faint with pn count 
rates of 0.5-7.9x10“^ ctss“^ and MOS count rates of 
0.1-2.6x10“^ ctss“^. Comparing their count rates in 
the three observations, we investigated their variabil¬ 
ity. Out of 17 sources, 16 of them were non-variable 
(< 3 (t). Only one source, located at R.A.=18^08™54"22, 
Decl.=—37°06'5074 (equinox J2000.0; Icr positional un¬ 
certainty is 074), exhibited variability in the count rate 
at a signihcance level of ^ 4 (t and ^ 3(7 between 2006 and 
2007 pn data and 2001 and 2007 MOS data, respectively. 
However, this source did not show a signihcant variation 
(^ 2cr) between the 2006 and 2007 MOS data. 

We further studied this variable source by htting its 
spectra with different models. The count rates were low 
in the 2001 and 2006 observations, and probably due to 
this reason, an absorbed power law can provide a good 
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fit to the spectra. When the column density Nh was set 
as a free parameter, unphysically large power-law indices 
of 3-6 were favored. If we fixed Nh = 1.3 x 10^^ cm“^, 
the Galactic value toward the source direction (Dickey 
& Lockman 1990), the indices were lowered to 2.9-3. 6 , 
and the obtained MOS (absorbed) fluxes were in a range 
of 1.7-4.1x10“^^ ergcm“^s“^ and the 2006 pn flux was 
l-4lo;5 X 10“^"^. However the 2007 pn spectrum, the most 
significantly detected (with a count rate of 7.91±0.71 x 
10 “^ cts s“^) among the sources and observations, can not 
be fit with a single model such as a power law (reduced 

> 2; 20 degrees of freedom). We searched the SIMBAD 
database and USNO Bl.O catalog, no radio or optical 
counterparts (generally down to 20 mag in R band) were 
found. Given the properties, we suggest that this source 
is either a low luminosity X-ray binary or a background 
galaxy. The power law index values are too high for an 
AGN (e.g., Ulrich et al. 1997), the largest class among 
Fermi LAT sources (The Fermi-LAT Gollaboration 2015). 

4. Discussion 

Carrying out maximum likelihood analysis of more than 
6 -year Fermi y-ray data of the source region of SAX 
J1808.4—3658, we have detected a y-ray source with 
the best-fit position consistent with that of AMXP. The 
source’s y-ray spectrum can be described by an exponen¬ 
tially cutoff power law. The obtained parameters of F = 
1.6±0.4 and Ec = 5.5±3.7 GeV are within the parameter 
ranges for pulsars (0.4 < F <2, 0.4 GeV < Ec < 5.9 GeV; 
see the Fermi second pulsar catalog, Abdo et al. 2013), 
although the uncertainties are large due to low counts of 
the source. In addition, a possible orbital modulation has 
also been detected. These results support the association 
of the y-ray source with SAX J1808.4—3658. 

Observational studies of the transitional MSP binaries 
(i.e., J1023+0038 and XSS J12270-4859) have shown that 
y-ray emission is brighter during their active state when 
an accretion disk appears (Stappers et al. 2014; Takata 
et al. 2014; Xing & Wang 2014) than that in the disk- 
free state, while the radio pulsars are possibly still ac¬ 
tive but not observable. Very likely in the latter state, 
the y-ray emission arises from the magnetosphere of the 
MSPs (e.g., Takata et al. 2014). In the former state, it 
has been suggested that either the y-ray emission is en¬ 
hanced due to inverse Compton (IC) scattering of a cold 
pulsar wind off the optical/infrared photons from the ac¬ 
cretion disk (Takata et al. 2014), or alternatively self¬ 
synchrotron Compton processes at the magnetospheric re¬ 
gion of a propellering neutron star is the possible working 
mechanism for producing y-ray emission (Papitto et al. 
2014). Considering the similarities between the quiescent 
state of SAX J1808.4—3658 and the active state of the 
transitional MSP binaries, it is likely that a same emis¬ 
sion mechanism also works in the AMXP and thus the 
observed y-ray emission is expected. 

The orbital modulation from this source was possibly 
detected. However because of the relative low significance 
and the unique modulation profile, we can not draw a 
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certain conclusion. The modulation has two brightness 
peaks around the inferior conjunction (phase 0.25, when 
the companion is in front of the neutron star) and the 
superior conjunction (phase 0.75, when the companion is 
behind the neutron star) respectively. Such modulation 
has not been observed in other MSP binaries. Usually 
there is only one brightness peak, either around the infe¬ 
rior conjunction (see, e.g. PSR J1023-I-0038 in Bogdanov 
et al. 2011 and XSS J12270-4859 in Xing & Wang 2014) 
and possibly due to the occultation of the photon emitting 
region by the companion, or around the superior conjunc¬ 
tion (see, e.g., PSR B1957-f20 in Wu et al. 2012, 2FGL 
J0523.3—2530 in Xing et al. 2014) and possibly due to the 
viewing angle of the intrabinary interaction region (Wu 
et al. 2012; Bednarek 2014). Moreover, the orbital signals 
in these MSP binaries are only seen when accretion disks 
are not present (Bogdanov et al. 2011; Bogdanov et al. 
2014; Xing & Wang 2014). Thus the possible orbital mod¬ 
ulation needs further confirmation from different studies. 
Phase-resolved spectra may help identify the spectral dif¬ 
ferences and the origin of the modulation. Unfortunately 
the photon counts from this y-ray source were too low to 
allow such analysis. 

Gonsidering that the y-ray emission is from SAX 
J1808.4—3658, the >0.1 GeV y-ray luminosity of the 
source is ^ 5.7(i§ 5 x 10^^ergs“^ (for the exponentially 
cutoff power-law model) at the source distance of 3.5 kpc 
(Galloway & Gumming 2006). The spin-down luminosity 
Esd of SAX J1808.4—3658 is ^ 9 x 10^^ ergs“^ (Hartman 
et al. 2009), indicating a y-ray conversion efficiency 
of 63%. The efficiency is above the ‘death line’ defined 
in Xing & Wang (2013) with the characteristic age of 
~ 12 X 10® yr for SAX J1808.4—3658 (calculated from the 
pulsar parameters given in Hartman et al. 2009), which 
supports the suggestion that older MSPs tend to have 
higher rj^ values. 

AGNs are the major class of Fermi LAT sources (The 
Fermi-LAT Gollaboration 2015), and they may be iden¬ 
tified from their strong variability (e.g., Ulrich et al. 
1997; Williamson et al. 2014). We have analyzed three 
sets of XMM-Newton X-ray imaging data of the SAX 
J1808.4—3658 field, and found 17 faint X-ray sources in 
the 2a error circle of the y-ray source. Among them, only 
one had 3-4cr low flux variations. However, this variable 
source did not have AGN-like emission, not supporting 
that it could be a background AGN. We caution that from 
the X-ray variability study, none of the other 16 sources 
are likely an AGN, but we can not totally exclude the 
possibility. In order to identify their nature from spectral 
properties, deep X-ray observations are needed. 

As we write the paper, the Fermi third source catalog 
is released (The Fermi-LAT Collaboration 2015), and we 
note that the source 3FGL J1808.4—3703 is reported to be 
detected at the region of SAX J1808.4—3658. The catalog 
position of 3FGL J1808.4-3703 is R.A. = 272P12, Deck 
= —37P05 (equinox J2000.0), consistent with the best-fit 
position we obtained within uncertainties (see Figure 1). 
Thus our data analysis is confirmed by the catalog results. 
The catalog source is also identified not to be unassociated 
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with any known type of objects in the recently available 
catalogs. 
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Table 1. Binned likelihood analysis results for SAX J1808.4—3658. 


Spectral model 

>0.2 GeV Flux 
(10“® photon cm“^ s“^) 

r 

Ec 

(GeV) 

TS 

Power law 

3.2 ± 0.9 

2.2 ± 0.1 


32 

Power law with cutoff 

2.4 ± 0.9 

1.6 ± 0.4 

5.5 ± 3.7 

37 
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Table 2. Flux measurements of SAX J1808.4—3658. 


E 

(GeV) 

Fiow/10-^" 

(erg cm“^ s“^) 

TS 

0.15 

1.1 

0 

0.33 

1.7 

0 

0.74 

1.4 ± 0.5 

10 

1.65 

1.0 ± 0.3 

10 

3.67 

0.8 ± 0.3 

9 

8.17 

0.5 ± 0.3 

5 

18.20 

0.5 ± 0.4 

4 

40.54 

0.6 

0 

90.27 

1.8 

0 

20.10 

5.6 

0 


Columns 2 and 3 list the energy 
flux (E2 X dN/dE) and the TS 
value in each energy bin, respec¬ 
tively. The fluxes without uncer¬ 
tainties are upper limits. 
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Fig. 1. TS map of a 2° x 2° region, with an image 
scale of 0.04° pixel" centered at the position of SAX 
J1808.4—3658 in >0.5 GeV energy range. The color 
bar indicates the TS values. The dark cross marks 
the optical position of SAX J1808.4—3658. The dark 
dashed and solid circles are the 2cr error circles of the 
best-fit position obtained by us and given in the Fermi 
third source catalog for 3FGL J1808.4—3703, respectively. 
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Fig. 2. 7 -ray spectrum of SAX J1808.4—3658. The 
exponentially cutoff power-law and the power-law fits 
obtained from maximum likelihood analysis are shown 
as the dashed curve and dot-dashed line, respectively. 
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Orbital Phase 

Fig. 3. 2-300 GeV light curve folded with using the X-ray 
orbital parameters (Hartman et al. 2009). The phase zero 
is at the ascending node of the MSP in SAX J1808.4—3658. 
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Fig. 4. XMM-Newton pn image of the SAX J1808.4—3658 
field. The large circle indicates the 2(T error circle 
of the Fermi 7 -ray source, in which 17 X-ray sources 
were detected (marked with small blue circles). The 
source found with 3-4cr flux variations is marked with V. 




